Thermal vs quantum decoherence in double well trapped Bose-Einstein condensates 
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The quantum and thermal fluctuations of the phase are 
investigated in a cold Bose gas confined by a double well 
trap. The coherence of the system is discussed in terms of 
the visibility of interference fringes in both momentum and 
coordinate space. The visibility is calculated at zero as well 
as at finite temperature. The thermal fluctuations are shown 
to affect significantly the transition from the coherent to the 
incoherent regime even at very low temperatures. The coher- 
ence of an array of multiple condensates is also discussed. 

After the first interference measurements carried out 
on two expanding and overlapping condensates jij the 
problem of the relative phase between Bose-Einstein con- 
densates has stimulated extensive theoretical work in 
the last few years (see for example and references 
therein) as well as new experimental attempts to point 
out Josephson like phenomena H|| . A challenging ques- 
tion is the transition from a coherent to an incoherent 
regime associated with the increase of the fluctuations of 
the relative phase between condensates confined in differ- 
ent traps. The transition can be in principle controlled 
by changing the tunneling probability between the wells 
either by tuning the height and the width of the bar- 
rier or the number of atoms. Most of theoretical works 
have focused so far on the ideal case of zero tempera- 
ture where the fluctuations of the phase are of quantum 
nature. These works have shown that if the tunneling 
probability is sufficiently small the ground state of an 
atomic gas interacting with repulsive forces is not glob- 
ally coherent, but exhibits new features of quantum na- 
ture, associated with the appearance of number squeezed 
configurations (see for example ||). 

The purpose of this letter is to discuss the relative 
importance of the quantum and thermal fluctuations in 
driving the transition between the coherent and incoher- 
ent regimes. Due to the low energy scale associated with 
the Josephson oscillation of the condensates the thermal 
effect may actually become crucially important even if 
one works at very low temperature. 

Let us consider a dilute gas of atoms confined by an 
external potential characterized by a double well along 
the x direction. The problem of an array of multiple 
wells will be discussed in the last part of the work. Near 
the bottom of the two traps the potential can be ap- 
proximated by a harmonic function characterized by the 



oscillator frequency ujho which provides the relevant fre- 
quency scale of the collective excitations of each conden- 
sate (actually, due to the 3D nature of the problem, the 
harmonic expansion will in general introduce three dif- 
ferent frequencies). We start our discussion by recalling 
the classical problem of the Josephson oscillation in the 
framework of the Gross-Pitaevskii (GP) theory for the 
order parameter. If the overlap between the condensates 
confined in the two wells is small we can naturally con- 
struct an approximate solution of the time dependent GP 
equation in the form 

* (r,t) = [*, (r, N{) e^' + * r (r, N r ) e^] , (1) 

where >Fz ir (r, Ni >r ) are real functions satisfying the sta- 
tionary GP equation in the left (I) and right (r) wells 
respectively and normalized to J dr"i>f r — Ni :T with 
Ni + N r = N. The time dependence of the solution ([!]) is 
contained in the phases <I>i jr and in the number of atoms 
Ni jT confined in each well. The equations of motion can 
be written in the canonical form d^/dt = dHj/d(Tik) 
and d(hk)/dt = -dH.j/d§ where 
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is the so called Josephson Hamiltonian |||7]] depending on 
the conjugate variables fik and $ where k = (Ni — N r ) /2 
is determined by the relative number of atoms in the two 
condensates and $ = <I>; — $ r is their relative phase. 
In deriving eq.(|^) we have assumed k <C N and in- 
troduced the relevant parameters E c — 2dfJ,i/dNi and 
Ej = (h 2 /2m) J dydz^^r/dx - ^ r d^i/dx] x=Q where 
Hi is the chemical potential of each condensate and all the 
quantities are evaluated at Ni — N r — N/2. The term in 
E c accounts for the interaction effects in each condensate 
and vanishes in the non interacting gas since in this case 
the chemical potential does not depend on the number of 
atoms. In the Thomas-Fermi limit instead the parame- 
ter E c takes the value E c = (4/5)/i;/iV ; due to the N 2 / 5 
dependence of the chemical potential of a Bose-Einstein 
condensate trapped by a 3D harmonic potential. In the 
Thomas-Fermi regime the value of Ec hence decreases 
with N. The term in Ej, which describes the tunneling 
probability between the two wells and is at the origin of 
the Josephson current, instead increases with N. Expres- 
sion (pi) for Hj can be easily generalized to include the 
effects of gravity. 



Hj = -E c k 2 



Ej cos $ . 



(2) 
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Equation ([j]) represents a good approximation to the 
solution of the GP equation only if the Josephson motion 
is decoupled from the other modes of the condensate. 
This implies that the quantities $ and k should vary 
slowly in time with respect to the typical time l/u>ho 
characterizing the internal motion of each condensate. 
The Hamiltonian (Q) corresponds to the classical problem 
of the pendulum. For small oscillations near equilibrium 
(k = 0, $ = 0) the motion is of harmonic nature and is 
characterized by the classical plasma frequency 



= - \J EjE c 
n 



(3) 



so that, in order to ensure the decoupling from the 
internal oscillations of the condensate, the inequality 
lu c i <C Who should be satisfied. 

Let us now proceed to quantize the classical Hamilto- 
nian (||). This is achieved by replacing the conjugated 
variables $ and Tik with operators satisfying the com- 
mutation relation [$, hk] = ih. It is convenient to work 
in the " ^-representation" where k — —id/d& and the 
quantum problem can be described by the Hamiltonian 



Hj = - 



Ecd*_ 

2 o<s> 2 



- Ej cos $ 



(4) 



acting in the space of the periodical functions of period 
2tt. The quantization introduces quantum fluctuations 
in the equilibrium state of the system whose nature de- 
pends on the ratio Ec/Ej. Because of the periodicity 
constraint the uncertainty relation obeyed by the fluc- 
tuations of $ ad A: can be conveniently written in the 
form 



(A/fc 2 )(Asin 2 $) > (cos$) 2 /4 



(5) 



where (AA 2 ) = (A 2 ) — (A) 2 is the variance of the general 
observable A. The uncertainty relation (||) reduces to 
the usual form (Afc 2 )(A$ 2 ) > 1/4 only if the phase is 
localized around <E> = Wt . The quantity 



(cos $) 



(6) 



characterizing the right hand side of the uncertainty in- 
equality (|B|) has a physical meaning of first importance. 
It provides the degree of coherence of the configuration 
and will be called the coherence factor. If the value 
of the phase is localized around zero, the value of a is 
close to unity (full coherence) and one can expand eq. @ 
to calculate the quadratic fluctuations of the phase as 
(A$ 2 ) = 2 (1 — a). If instead the phase is delocalized 
and all its values are equally probable, then the value of 
a is zero (absence of coherence) . 

The parameter a is directly related to measurable 
observables. For example if one calculates the in situ 
momentum distribution of two separated condensates 
with a fixed relative phase <E>, one finds fringes in mo- 
mentum space given by the expression n (p) = 



2 [1 + cos (p x d/h + $)] rio(p) where no(p) is the momen- 
tum distribution of each condensate and d is the relative 
distance along the x axis. The typical width of no(p) 
is given by h/R where R is the spatial size of each con- 
densate. The momentum distribution of a single trapped 
Bose gas has been already the object of measurements 
based on stimulated light scattering |l(|. If the value of 
the relative phase fluctuates the average value of the mo- 
mentum distribution at equilibrium is described by the 
equation 



n (p)) = 2 [1 + a cos (p x d/K)] n ; (p) 



(7) 



where the parameter a corresponds to the average (cos $) 
and hence coincides with the coherence factor introduced 
in eq. (|J) (we have used the fact that (sin $) = at equi- 
librium). One should not however confuse the parameter 
a with the visibility of fringes obtained in a single real- 
ization of the experiment which can give rise to a well 
defined interference pattern even if the initial state is 
not coherent. The parameter a corresponds instead to 
an average taken on several measurements and is conse- 
quently an intrinsic property of the state of the system. 
This distinction is crucial if one considers interference 
experiments with two condensates. A similar structure 
characterizes also the fringes of the density distribution 
measured after expanding the condensates |Q. In this 
case in order to obtain a simple result one has to make 
the additional assumption that the condensates do not 
interact each other when they overlap. For long expan- 
sion times t the phases of the two condensates evolve 
according to the law (r±d/2) 2 m/M jll| so that the full 
density profile takes the following asymptotic form 



n(r,t)) 



2a cos 



nt 



(8) 



where n± = no(r ± d/2,t) are the densities of the two 
expanding condensates. Eqs.(Qj8|) show that the value 
of a can be inferred from either the measurement of the 
momentum distribution and/or of the density profile af- 
ter expansion, although in the second case the conditions 
of applicability of eq.(^) are more severe. 

Let us now discuss the behaviour of the ground state 
of the Josephson Hamiltonian (^) . In the limit of strong 
tunneling Eq j Ej <C 1 the system undergoes small os- 
cillations around the equilibrium value $ = and the 
Hamiltonian becomes quadratic in <F In this limit 
the fluctuations of the phase are given by (A$ 2 ) = 
(y/Ec/Ej)/2 <C 1 and are hence small. One can con- 
sequently regard the inequality Eq /Ej <C 1 as the phys- 
ical condition for applying the classical GP theory to the 
problem of the double-well potential and, consequently, 
to treat the systems as a globally coherent object de- 
scribed by a unique order parameter. Also the fluctua- 
tions of k in the ground state can be easily calculated 
and one finds the result (Afc 2 ) = ( v / E J /E c )/2 > 1. In 
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the strong tunneling limit the uncertainty relation (||) re- 
duces to the usual form (Afc 2 )(A$ 2 ) > 1/4 and actually 
becomes an identity as a consequence of the harmonic na- 
ture of the Hamiltonian. Notice, however, that if (A$ 2 ) 
is of the order of 1/N or smaller, the formalism devel- 
oped above is no longer adequate and should be improved 
through both the inclusion of the factor yl — 4fc 2 /A 2 in 
the term — Ej cos <$> of the Josephson Hamiltonian (Q), 
and a more microscopic approach to the phase operator 
(see for example and references therein). In the fol- 
lowing we will mainly limit our discussion to the case 
Ec S> Ej/N 2 which ensures that the quadratic fluctua- 
tions of the phase are in any case larger than 1/N. 

In the case of weak tunneling Ec/Ej ^> 1 the be- 
haviour of the fluctuations is very different. The Joseph- 
son term Ej entering the Hamiltonian (^) can be ne- 
glected in first approximation and the wave function of 
the ground state is simply a constant given by 1/V2n, 
showing that the relative phases between the two con- 
densates is distributed in a random way. At the same 
time the fluctuation of the relative number k of atoms 
in the two traps becomes smaller and smaller and vanish 
like 2(Ej / Ec) 2 ■ This is consistent with the uncertainty 
relation (||) since in the same limit the coherence fac- 
tor vanishes like a — ► 2(Ej/Ec)- In the weak tunneling 
limit the ground state of the system exhibits fragmented 
Bose-Einstein condensation and is characterized by the 
macroscopic occupation of the two single-particle states 
localized in the two wells. 

In fig.l we show the coherence factor a as a func- 
tion of the ratio Ec/Ej calculated by solving explicitly 
the Schrodinger equation for the ground state with the 
Josephson Hamiltonian (Q). The figure shows that for 
values of Ej smaller than Ec the coherence factor is sig- 
nificantly quenched, pointing out the occurrence of a con- 
tinuous transition to the number-squeezed regime. The 
transition is accompanied by a change of the lowest exci- 
tation energy ui p of H j from the classical value (|J) to the 
"free" value Ec/2 which is obtained by setting Ej = 
in (Q) . An accurate expression for u> p is obtained by eval- 
uating the ratio between the cubic and energy weighted 
sum rules ([[k, Hj],[Hj,[Hj,k]]}) = EjE c (cos<l> 2 } and 
([k , [Hj, k))) = Ej(cos $), relative to the operator k. The 
result is 



1 / (cos 2 g) 

h\ (cos <P) 



(9) 



where the averages should be evaluated on the ground 
state. In the limit Ec <C Ej one has (cos 2 $) = (cos <I>) = 
1 and one recovers the classical result (^). In the opposite 
limit one finds (cos$ 2 ) = 1/2, (cos$) = 2Ej/E c and 



E c /2h. 



So far we have investigated the problem in an ideal 
situation at zero temperature. Due to the smallness 
of the plasma frequency (p|) it is important to under- 



stand the role of the thermal fluctuations of the rel- 
ative phase of the two-condensate system. These are 
expected to become important as soon as the temper- 
ature is of the order or higher than Tiuj p . To investi- 
gate the thermal effect we have calculated the thermal 
average a(T) = £ n a n exp[-£„/T]/ exp[-E n /T\ of 
the coherence factor where n and E n are the eigen- 
states and eigenenergies of the Josephson Hamiltonian 
(0), a n = (n | cos$ | n) are the corresponding quan- 
tum averages and we have set the Boltzmann constant 
equal to 1. The results are presented in fig. 2 as a func- 
tion of the parameter T/Ej for two different values of 
Ec / Ej . At T = one recovers the values of a given in 
fig.l. One clearly see that even if quantum effects are 
small the thermal decoherence of the phase becomes im- 
portant at temperatures of the order T ~ Ej. The full 
line in the same figure gives the classical prediction 



a c i{T)-- 



d$ cos $ exp [Ej cos $/T] 
J^d$exp[Ej cos$/T] 



(10) 



holding for T > cu p . If T < Ej, eq. (|10|) gives the result 
(A$ 2 ) = T/Ej <C 1 for the classical fluctuation of the 
phase. The fluctuation of k can be also easily calculated 
starting from the classical expression @ and one gets the 
result (Afc 2 ) = T/Ec for any value of T. The thermal 
fluctuations of k may be large and consequently at finite 
T the system will not in general exhibit the phenomenon 
of number squeezing, even if the condition Ej <C Ec is 
satisfied. 

It is useful to discuss the relevant scale of energies and 
temperatures in a specific example. Let us consider a 
system of 10 3 atoms per well and the value fi = bfiLUho 
for the chemical potential. The coupling constant E c 
is about 5 x 10~ 3 hu>ho- By choosing Ej = Tiuiho one has 
lu p ~ O.lujho- With the above choices both the conditions 
lu p <C ujho and Ec Ej/N 2 needed to apply the for- 
malism developed above are well satisfied. At T = the 
system is practically coherent (a ~ 1) since Ec -C Ej. 
However if T is of the order of the oscillator energy hojho, 
which corresponds to a rather low temperature in stan- 
dard experiments, the coherence is partially lost because 
of thermal effects. By changing the value of Ej and/or 
T one can clearly explore other interesting regimes. 

The Josephson effect with two condensates can be nat- 
urally extended to an array of multiple wells. This ex- 
tension is particular relevant because of the recent ex- 
perimental efforts to investigate one-dimensional optical 
lattices @. The array of condensates confined in multi- 
ple wells is described by the multiple-condensate Hamil- 
tonian 

£j = -jEc E^- £ .'E cos (^+i - ( n ) 



fc=l 



k=l 



where N s is the number of wells, is the phase of the 
k-th condensate and Ec, Ej are the parameters charac- 
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terizing the double- well Hamiltonian (0) . In the large N s 
limit the one-dimensional Hamiltonian ([ll]) is known to 
exhibit a phase transition at zero temperature occurring 
at the critical value E c = 1.62£ , , / Q. The "superfluid" 
phase is not accompanied by the appearence of an or- 
der parameter and cannot be described in terms of mean 
field theories. It is instead characterized by an algebraic 
decay of the phase correlation function (cos($fe — <!>;)) 
at large distances | k — I |^> 1 (quasi long range order). 
The decay becomes exponential for Ec > l.Q2Ej. The 
phase transition disappears at finite temperature. Actu- 
ally in the "classical" regime of high temperatures the 
phase correlation functions can be evaluated analytically 
and one finds the simple result 



(cos($ fe - $,)>c; = 4; 



(12) 



which exhibits an exponential decay for large values of 
| k — I \. In eq.([l2|) a c i is the classical value (1C) cal- 
culated for the double well problem. Result (12) can 
be used to evaluate the momentum distribution n(p) 
of the array. By taking the Fourier transform \&(p) = 
(2irh)~ 3 / 2 f c?r^E'(r) exp[— ip ■ r/h] of the order parame- 



ter \ff(r) 



4'fe(r)e k and using the relationship 



*fc(p) = ^i(p) exp[— i(fc — j)p x d/h] following from the 
translational invariance properties of the array, one finds 
the result 



n(p) =n (p)^e 



— i{k — l)p x d/U 



k.l 



N s n {p) 



1 



1 



cos($ fc -$,)> (13) 



(14) 



2a ci cos(p x d/h) 



where no(p) is the momentum distribution of each con- 
densate, d is the distance between two consecutive con- 
densates and, in deriving the last identity, we have as- 
summed N s » 1(1 — a c i). Result ( |l4| ) explicitly points 
out the effects of coherence. When cos(p x d/K) = 1 
the incoherent signal N s no(p) is amplified by the fac- 
tor (l + ay)/(l — a c i). Viceversa, when cos(p x d/h) = —1 
the signal is suppressed by the factor (1 — a c i)/ (1 + u c i)- 
In conclusion we have investigated the consequences 
of the quantum and thermal fluctuations of the phase 
on the coherence phenomena exhibited by Josephson like 
configurations. Thermal effects are predicted to play a 
major role even at very low temperature and should be 
consequently taken into account in order to control the 
transition to the new quantum phases exhibited by these 
ultracold Bose gases. 
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Fig.l Coherence factor a at zero temperature as a func- 
tion of the ratio Ec/Ej. 

Fig. 2 Coherence factor a as a function of the ratio 
T/Ej for E c /Ej = 3 (a), E c /Ej = 1 (b). The classical 
curve ( pF| ) is also shown (c). 
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